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We introduce a four-parameter extended family of distributions related to the wrapped Cauchy 
distribution on the circle. The proposed family can be derived by altering the settings of a 
problem in Brownian motion which generates the wrapped Cauchy. The densities of this family 
have a closed form and can be symmetric or asymmetric depending on the choice of the param- 
eters. Trigonometric moments are available, and they are shown to have a simple form. Further 
tractable properties of the model are obtained, many by utilizing the trigonometric moments. 
Other topics related to the model, including alternative derivations and Mobius transformation, 
are considered. Discussion of the symmetric submodels is given. Finally, generalization to a 
family of distributions on the sphere is briefly made. 
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1. Introduction 

As a unimodal and symmetric model on the circle, the wrapped Cauchy or circular 
Cauchy distribution has played an important role in directional statistics. It has density 

/(g;/i,p) = ^- 1 ~ P -, -7i<e<n, (1) 

2n 1 + p z — 2pcos(0 — (j,) 

where p (g [— 7t, 7t)) is a location parameter, and p (g [0,1)) controls the concentration 
of the model. This distribution has some desirable mathematical features, as discussed 
in Kent and Tyler [11] and McCullagh [15]. We write 9 ~ WC(/i, p) if a random variable 
9 has density (1). 
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In modeling symmetric circular data, the wrapped Cauchy distribution could be one 
choice as well as other familiar symmetric models such as the von Mises and wrapped 
normal distributions; see also Jones and Pewsey [7]. In reality, however, it is not so 
common to find such data in fields of application. In many cases the data of interest are 
asymmetrically distributed, and therefore probability distributions having asymmetric 
densities are desired. 

Construction of a tractable circular model with an asymmetric shape has been a prob- 
lem in statistics of circular data. To tackle this problem, some asymmetric extensions of 
well-known circular models have been proposed in the literature. Maksimov [12], Yfantis 
and Borgman [20] and Gatto and Jammalamadaka [4] discussed an extension of the von 
Mises distribution, generated through maximization of Shannon's entropy with restric- 
tions on certain trigonometric moments. Batschclct [2] proposed a mathematical method 
of skewing circular distributions that has seen renewed interest very recently. Pewsey [16] 
presented a four-parameter family of distributions on the circle by wrapping the stable 
distribution, which is asymmetric, onto the circle. Recent work by Kato and Jones [9] pro- 
posed a family of distributions arising from the Mdbius transformation which includes the 
von Mises and wrapped Cauchy distributions. Unlike familiar symmetric distributions, it 
is often difficult to deal with skew models in statistical analysis. This difficulty is partly 
due to the lack of some mathematical properties that many of the well-known symmetric 
models have. For example, existing asymmetric models often have complex normalizing 
constants and trigonometric moments, which could cause trouble in analysis. 

In this paper, we provide a four-parameter extended family of circular distributions, 
based on the wrapped Cauchy distribution, by applying Brownian motion which, to our 
knowledge, has not been used to propose a skew distribution on the circle. The four 
parameters enable the model to describe not only symmetric shapes but also asymmetric 
ones. An advantage of the proposed model is its mathematical tractability. For instance, 
it has a simple normalizing constant and trigonometric moments. The current proposal 
is more tractable than the family discussed by Kato and Jones [9] , but is complementary 
to it in the sense that the latter has other advantages (particularly some associated with 
Mdbius transformation). 

The subsequent sections are organized as follows. In Section 2 we make the main 
proposal of this paper. The derivation of the proposed model is given, and the probabil- 
ity density function and probabilities of intervals under the density are discussed. Also, 
special cases of the model are briefly considered. Section 3 concerns the shapes of the den- 
sity. Conditions for symmetry and for unimodality are explored, and the interpretation 
of the parameters is discussed through pictures of the density. In Section 4 we discuss the 
trigonometric moments and problems related to them. It is shown that the trigonometric 
moments can be expressed in a simple form. The mean direction, mean resultant length 
and skewness of the model are also considered. Some other topics concerning our model 
are provided in Section 5. Apart from the derivation given in Section 2, there are other 
methods to generate the family, which are discussed in that section. In addition, we study 
the conformal transformation properties of the distribution. In Section 6 we investigate 
some properties of the symmetric cases of the proposed model. These submodels have 
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some properties which the general family does not have. Finally, in Section 7, general- 
ization to a family of distributions on the sphere is made, and its properties are briefly 
discussed. 

Further material on the topic of this paper can be found in its technical report version, 
Kato and Jones [10]. This includes the path of the expected exit points of a Brownian 
particle, random variate generation using the Markov chain Monte Carlo method, and 
proofs of Theorems 1, 3 and 10. 

2. A family of distributions on the circle 

2.1. Definition 

It is a well-known fact that the wrapped Cauchy distribution can be generated as the 
distribution of the position of a Brownian particle which exits the unit disc in the two- 
dimensional plane; see, for example, Section 1.10 in Durrett [3]. In addition to the 
wrapped Cauchy circular model, some distributions on certain manifolds are derived 
from, or have relationships with, Brownian motion. The Cauchy distribution on the real 
line is the distribution of the position where a Brownian particle exits the upper half 
plane. Considering the hitting time of the particle, one can obtain the inverse Gaussian 
distribution on the positive half-line. Kato [8] proposed a distribution for a pair of unit 
vectors by recording points where a Brownian particle hits circles with different radii. 

By applying Brownian motion, we provide a family of asymmetric distributions on the 
circle, which includes the wrapped Cauchy distribution as a special case. The proposed 
model is defined as follows. 

Definition 1. Let {B t ;t > 0} be R 2 -valued Brownian motion without drift starting at 
Bo = Pi(cos/ii,sin/xi)' ; where < p\ < 1 and — n< /ii < n. This Brownian particle will 
eventually hit the unit circle. Let n be the first time at which the particle exits the 
circle, that is, t\ = inf{f; \\B t \\ = 1}. After leaving the unit circle, the particle will hit a 
circle with radius p^ 1 (0 < p2 < 1) first at the time T2, meaning Ti = inf{t; \\B t \\ = p^ 1 } ■ 
Then the proposed model is defined by the conditional distribution of B Tl given B T2 = 
P2 1 (cos H2, sin P2)' , where —Tt<p2 < 71. 

To put it another way, the proposed random vector, B Tl given B T2 , represents the 
position where a Brownian particle first hits the unit circle, given the future point at 
which the particle exits a circle with a larger radius. From the next subsection, we 
investigate some properties of the proposed model. 

2.2. Probability density function 

One feature of the proposed model is that it has a closed form of density with simple 
normalizing constant. It is given in the following theorem. See Appendix A of Kato and 
Jones [10] for the proof. 
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Theorem 1. Let be a [— 7t,7t) -valued random variable defined as B Tl = (cosO,sin0)'. 
Then the conditional density of given B n = p^ 1 (cos p 2 , sin p 2 )' is given by 

f(6) = Cl{l + pi-2p 1 cos(9- f i 1 )}{l + p 2 2 ~2p 2 cos(6~p 2 )}}-\ -n<9<n, (2) 

where —tc < p±, p 2 < 7t, < pi,p2 < 1> and the normalizing constant C is 

(1 - ~ P2H 1 + P1P2 ~ 2 PiP2 cos(^i - p 2 )} 



C = 



2ti(1 -pip*) 



For convenience, write ~ EWC(pi, p 2 , Pi, P2) if a random variable has density (2) 
( "E" for "Extended" ) . Note that the density does not involve any infinite sums or special 
functions. 

It is easy to see that distribution (2) reduces to the wrapped Cauchy (1) when either 
Pi or p 2 is equal to zero. 

Quite often, it is advantageous to express the random variable and parameters of the 
proposed model in terms of complex numbers rather than real numbers. Define a random 
vector by Z = e' e where ~ EWC(/ii, p 2l pi, p 2 ). Then Z has density 

J.H-^i-NM-l^ 

JK ' 27Tl-|0 1 2 |2|z-</> 1 |2|z-0 2 |2' ^ 

with respect to arc length on the circle, where <fii — pie lftl , <fi 2 — p 2 c lfX2 and 3D = {z£ 
C; \ z\ = 1}. It is clear that the parameters in this formulation, <f>i and (j> 2 , take values 
on the unit disc in the complex plane denoted by D = {z G C; |z| < 1}. For brevity, we 
denote the distribution (3) by EC*(^i, <fo). Also, as in McCullagh [15], write Z - C*(</>i) 
if Z follows distribution (3) with <j) 2 = 0. 

We remark that the above change-of-variable and reparametrization do not actually 
change the distribution. In later sections, we utilize either of the representations (2) or 
(3), whichever is the more convenient. 

2.3. Probabilities 

As well as the density of the proposed model, the probabilities of intervals under the 
density can also be expressed without using infinite sums or special functions. 

Theorem 2. Let a random variable follow the EWC(/Xi, p 2l pi, p 2 ) distribution. If 
Pi =/= P2 or p\ p 2 , then the probability of intervals under the density of is given by 

P(a<0<6) 

= [ c[{i + pi-2 Pl cos(e-p 1 )}{i + p 2 2 -2p 2 cos(e-p 2 )}}- 1 de 
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, 2/0i{pi(l + pi) - p 2 (l + Pi)cos(m - p 2 )} 



1-P? 

1 + pi (0- Pi 
arctan< tan 



i 6 



1 - Pi V 2 
2p 2 {p 2 (l + Pi) - Pi(l + pi) cos(^i - pt 2 )} 



i-pI 



M2 



l-/7 2 — \, 2 

where C is defined as in Theorem 1, —n < a < b < n, 

D = {p\+ pl)(l + p 2 p 2 ) - 2pip 2 (l + pf )(1 + p 2 2 ) cos(mi - M2) + 4p?p^ cos 2 (/ii - // 2 ) 



tan{i(a — p)} > tan{|(6 - /i)}, 
otherwise. 



If pi = p 2 and pi = pi, then 



P(a<@<b)= ( C{l + p\ -2 / 9iCos(6»-/ii)}- 2 d6» 

J a 



2C 



(i -pi) 2 



pi $m{U - pi, 



1 + p\ — 2pi cos(# — pi 
, 1 + 



arctan^tanr-^ 



I-Pi 
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Proof. The result is straightforward from equations (2.553.3), (2.554.3) and (2.559.2) of 
Gradshtcyn and Ryzhik [5]. □ 



2.4. Special cases 

The proposed family (2) contains some known distributions as special cases. 

Case 1: The wrapped Cauchy distribution. As mentioned in Section 2.2, model (2) 
becomes the wrapped Cauchy WC(^i,pi) if p 2 = 0. Similarly, the model is WC(p 2 , p 2 ) 
when pi = 0. 

Case 2: A special case of the Jones and Pewsey [7] family. If p\ = p 2 (= p) and p»\ = 
p 2 (= p), then density (2) reduces to 



M = 



i 



i i-/) 2 

2n 1 + p 2 \ 1 + p 2 - 2pcos(6 - p) 



—n < 9 < 7t. 



(4) 
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This submodel corresponds to a special case of the family presented by Jones and 
Pewsey [7]. Their model has density 

{cosh(^) + sinh(^)cos(6l-^)} 1 /^ „^ a ^„ 

/JP W = 7T~5 7 r 7 7yi > -7T<6*<7T, 

27tP 1 /^,(cosh(K?/!)) 

where k > 0, ^eR, — 7T < /i < 7t and Pi/^,(z) is the associated Legendre function of 
the first kind of degree 1/ip and order (Gradshtcyn and Ryzhik [5], Sections 8.7, 8.8). 
Our submodel is equivalent to their family with ip = — 1/2 and k = 21og((l + p)/(l — 
p)). In this case the associated Legendre function simplifies to P_2(z) = P\(z) = z. (Sec 
equations (8.2.1) and (8.4.3) of Abramowitz and Stegun [1].) This fairly heavy-tailed 
circular distribution is shown by Jones and Pewsey to tend, suitably normalized, to the 
t distribution on V3 degrees of freedom as (/? — > 1 and hence) k — > oo! 

Case 3: One-point distribution. As p\ {p2) tends to one, the model converges to a 
point distribution with singularity at 9 — p± (6 = p2). Normalizing by a scale factor of 
(1 — p)/yfp, the limiting version of the wrapped Cauchy distribution as p — > 1 is the 
ordinary Cauchy distribution. That argument can be extended to show that the Cauchy 
limiting distribution also arises in this case. 

Case 4-' Two-point distribution. Assume that p\ = p-2 (= p) and piy^p-2- When p goes 
to one, model (2) converges to the distribution of a random variable which takes values 
on pi or P2 with probability 0.5. If p\ < P2 arid p tends to one, {^/p/{l — p)}(Q — 
pi) converges to ^C(0, 1) + \l, where C(0, 1) is the standard Cauchy, and / has the 
distribution function F(x) = limA-i-oo l/[7t{l + (f — A) 2 }] dt. Roughly speaking, this 
limiting distribution is a 50:50 mixture of the standard Cauchy and a point distribution 
with singularity at X = 00. Similarly, one can show that if p\ > P2 and p tends to one, 
{y/p/(l — p)}(© — Pi) converges to a 50:50 mixture of the standard Cauchy and a point 
distribution with singularity at X = —00. The limiting distribution of {^fp/ (1 — p)}(6 — 
PL2) can be discussed in a similar manner. 

Case 5: The circular uniform distribution. When p\ = P2 = 0, the distribution reduces 
to the circular uniform distribution. 



3. Shapes of probability density function 
3.1. Conditions for symmetry 

As briefly stated in Section 1, the proposed model can be symmetric or asymmetric, 
depending on the values of the parameters. The condition for symmetry is clearly written 
out in the following theorem. 

Theorem 3. Density (2) is symmetric if and only if p\ = P2, Pj = or p\ = P2 + (j — 1)tc 
(j = l,2). 

See Appendix B of Kato and Jones [10] for the proof. The three-parameter symmetric 
submodels arising under the above conditions will be discussed in Section 6. 
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3.2. Conditions for unimodality 

Turning to conditions for unimodality, it is possible to express these in terms of an 
inequality. The process to obtain the inequality is similar to that in Kato and Jones [9], 
Section 2.5. In the following discussion, take p 2 = without loss of generality. First we 
calculate the first derivative of density (2) with respect to 9, which is given by 

d [pi sin(fl -m)(l + pl- 2p 2 cosfl) + p 2 sinfljl + p\ - 2p x cos(6> - p x )}\ 

d9 I[ jOC {l + pj-2p 1 cos(9-p 1 )} 2 (l + p 2 -2p 2 cos9) 2 

Then it follows that the extrema of the density are obtainable as solutions of the following 
equation: 

ao + ai cos 9 + a 2 sin 9 + 0,3 cos 9 sin 9 + 04 cos 2 9 = 0, (5) 

where 

a = 2/>i/9 2 sin / ui, a± = p±(l + p 2 )smp 1 , a 2 = -pi(l + pi) cos^i - p 2 (l + p\), 
a 3 = 4pi/92 cos //i, a 4 = — Apxp 2 sin/Lti. 

Following Yfantis and Borgman [20] and Kato and Jones [9], Section 2.5, put x = tan(0/2) 
so that cos6* = (1 — x 2 )/(l + x 2 ), s'm9 = 2x/(l + x 2 ). It follows that equation (5) can 
be expressed as a quartic equation in x whose discriminant D (Uspensky [19]) can be 
written down in terms of clq, . . . , 04 as in (6) of Kato and Jones [9]. The quartic equation 
has four real roots or four complex ones if D > 0, and two real roots and two complex 
ones if D < 0. Therefore the distribution is bimodal when D > and unimodal when 
D < 0. Since aj (0 < j < 4) are functions of p,\,p\ and p 2 , we can write the conditions 
for unimodality as a function of these three parameters. For general p 2 , it is easy to see 
that the discriminant is expressed in terms of p\ — pi, p\ and p 2 . 

Figure 1 exhibits areas of positivity (bimodality of density) and negativity (unimodality 
of density) of the discriminant D for two pairs of parameters. Figure 1(a) suggests that 
density (2) with p\ = 0.5 becomes unimodal for most cases when —tc/2 < pi < rc/2. It 
seems that bimodality is most likely to occur if tc/2 < pi < 3tc/2. The other frame of 
Figure 1 seems to show that distribution (2) with p 2 = 2tc/3 is bimodal if p\ or p 2 
is sufficiently large, and p 2 j p\ takes a value close to 1. In particular, when p\ = p 2 , 
the range of parameters corresponding to unimodality seems less wide than that when 
pi = ap 2 with a 7^ 1. As will be shown in Section 6.2, for p\ = 27t/3 and p\ = p 2l the 
density becomes unimodal if p\ < 2 — a/3 — 0.268 and bimodal otherwise. 

3.3. Pictures of density 

We plot density (2) for selected values of the parameters in Figure 2. Confirming the 
results in the previous two subsections, this figure shows that the density can be sym- 
metric or asymmetric, unimodal or bimodal, depending on the choice of the parameters. 
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-1.0 -0.5 0.0 0.5 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



rho1*cos{mu1) rhol 

(a) (b) 

Figure 1. Discriminant D as functions of: (a) (pi cos //i, pi sinpi) with p2 = 0.5, and (b) (pi,p2) 
with pi = 27t/3. Positive D, corresponding to bimodality, is displayed in black, and negative D, 
corresponding to unimodality, is exhibited in white. 

Figure 2(a), displaying the density for fixed values of p,2,pi and P2, suggests that the 
density is symmetric if pi = P2 or pi = p% + n and asymmetric otherwise. Figure 2(b) 
exhibits how the shape of density (2) changes as p\ increases. It is clear from this frame 
that the greater the value of p\ , the greater the concentration of the density. The frame 
also implies that the absolute value of the circular skewness of the density with fixed \i\ 
and H2 takes a large value if P1/P2 is close to 1. Some symmetric cases of the density 
are shown in Figure 2(c). It appears to be that the density is unimodal when /Lti = p2i 
whereas bimodality can occur when (12 = P>x + 7t. Figure 2(d) is a case in which two 
parameters, p\ and p2, arc fixed to be equal (to 0.5), also implying that density (2) is 
symmetric. We will focus on this submodel in Section 6.2. 

4. Trigonometric moments and related problems 
4.1. Trigonometric moments 

It is often the case that asymmetric distributions on the circle have complicated trigono- 
metric moments or Fourier coefficients. One feature of our asymmetric model, however, 
is the relative simplicity of its trigonometric moments. The expression for the moments 
is greatly simplified if the variables and parameters are represented in terms of complex 
numbers. 
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Theta Theta 
(c) (d) 

Figure 2. Density (2) for p2 = and (a) pi = 2/3, p2 = 1/3 and: pi = (solid), 7t/3 (dashed), 
27t/3 (dot-dashed) and 7t (dotted), (b) pi = rc/2, p 2 = 1/3 and: pi =0 (solid), 1/4 (dashed), 1/2 
(dot-dashed) and 3/4 (dotted), (c) p 2 = l/3 and: (pi, pi) = (7t, 1/3) (solid), (0,0) (dashed), 
(0,1/3) (dot-dashed) and (0,2/3) (dotted), and (d) pi = p 2 = 0.5 and: pi = (solid), tt/3 
(dashed), 27t/3 (dot-dashed) and n (dotted). 



Theorem 4. Let a random variable Z have the EC*(<^i, </>2) distribution. Then the nth 
trigonometric moment of Z is given by 



E(Z n ) = 



(1 - |0 2 | 2 )(1 - - (1 -^i | 2 )(1 - 0i<fe)^ +1 , 

(0i-^ 2 )(l-|0i^| 2 ) ' 

l + n + (l_ n )|0 1 |2 

TTRTp 01 ' ^ 1 = 02 ' 
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Proof. The nth trigonometric moment can be expressed as 



E(Z") = 




If <fii 7^ 02, the integrand in (6) is holomorphic on the unit disc, except at two poles of 
order 1, that is, z = 0j, j = 1,2. When 0i = 02, the integrand in (6) has a single pole 
of order 2 at z = <f>\. Hence, from the residue theorem (e.g., Rudin [18], Theorem 10.42), 
we obtain Theorem 4. □ 



4.2. Mean direction and mean resultant length 

Let Z be distributed as EC*^!,^)- From Theorem 4 it is easy to see that the first 
trigonometric moment for Z is simply expressed as 

w(7\ a-i<fei 2 )^i + (i-i0ii 2 )<fe f n , (7 , 

E(Z) = for any 0i and 02- (7) 

1 - 1 <Pi 02 r 

The following result provides the condition under which the first trigonometric moment 
is equal to zero. 

Corollary 1. The necessary and sufficient condition for E{Z) = is given by 0i = —02- 

Proof. It is straightforward to show the sufficient condition for E(Z) = 0. The necessary 
condition is proved as follows. Let E(Z) = and 0i ^ —02- Then, from Theorem 4, 
the following equation holds between the two parameters: 0102 = (0i +02)/(0i + 02)- 
Taking the absolute values of both sides of the above equation, we have |0i02| = 1- Since 
\4>j\ < 1 (j = 1,2), there does not exist 0i and 02 such that |0i02| = 1- Therefore, if 
0i ^-02, then £(Z) ^0. □ 

The mean direction is a measure of location which is defined as &rg{E(Z)} if E(Z) ^ 
and is undefined if E(Z) = 0. The mean resultant length, which is a measure of concentra- 
tion, is defined by |i?(Z)|. Figure 3 shows the mean direction and mean resultant length 
of a random variable having density (3) with fixed 02 = 0.5. Figure 3(a) suggests that 
the mean direction monotonically increases as arg(0i) increases from —n to n. As seen in 
Corollary 1, Figure 3(b) confirms that the mean resultant length is zero if 0i = — 0.5. It 
seems from this frame that unimodality holds for the mean resultant length as a function 

of 01. 
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Figure 3. Plot of (a) mean direction and (b) mean resultant length of a random variable having 
the distribution EC*(0i,O.5), as a function of (f>i. 



4.3. Skewness 

A measure of skewness for circular distributions (Mardia [13]), s, is denned by 

s = E[lm{{Zc~^) 2 }]/{l-8f' 2 , (8) 

where £ and 5 are the mean direction and mean resultant length for a 9-D-valucd random 
variable Z , respectively. For our model, it is possible to express the skewness for the 
proposed model in a fairly simple form as follows. 

Corollary 2. Let Z be distributed as EC*(</>i, fa). Then the skewness for the distribution 
of Z is given by 



= 11 rr=fc >- 2 (i - P )- 3/2 Im(fafa)(\^\ 2 - \h\ 2 )(i - I0il 2 )(i - \h\ 2 ), 
(i - mfar 



where p=\(l- \fa\ 2 )fa + (1 - \fa\ 2 )fa\/(l - \fafa\ 2 ). 
Proof. Note that the numerator of (8) can be expressed as 



E[lm{(Zc~^r}]=Im[E(Z 2 ){E(Z)/\E(Z)\n 

Then, by using Theorem 4 and equation (7), the skewness (9) can be obtained after 
lengthy, but straightforward, calculations. □ 
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Now we write s = s(0i,02) m order to view the skewness (9) as a function of 0i and 
02. It follows from Corollary 2 and Theorem 3 that the following properties hold for s. 

Theorem 5. The following properties hold for the skewness (9): 

1. s(0i, <p2) = •<=>■ density (3) is symmetric. 

2. s(0i,^) >0 (<O)j^Im(0i^)(|^i| - 1021) >0 (<0). 

3. s(01,0 2 ) = -s(01,0 2 0?/|0l| 2 ) = - S (010|/|02| 2 ,02),01,02^O. 

4. s(02,0i) = s(0i,0 2 ). 

5. s(0i,0 2 ) = -s(0i,0 2 ). 

6. s(a0i, a0 2 ) = s(0i, 02), a <E <9ZX 

7. Assume p\,p% > and sin(/^i — /X2) > (< 0). Then lim pi ^i s(pie 1AI1 , P2C 11 * 2 ) = 
00 (—00). 

Figure 4 plots this skewness when 02 = 0.5. As the first property of Theorem 5 shows, 
the skewness is equal to zero if and only if Im(0i0 2 ) =0, |0i| = |02| or \c/)j \ = (j = 1, 2). 
This figure also confirms the second property of Theorem 5 that the model is positively 
(negatively) skewed if |0i| > |02| and Im(0i0 2 ) > (Im(0i0 2 ) < 0), or |0i| < |0 2 | and 
Im(0i02) < (Im(0i02) > 0). It can also be seen that the fifth and seventh properties of 
Theorem 5 hold in this figure. 




-1.0 -0.5 0.0 0.5 1.0 



Re(phi1) 

Figure 4. Skewness (9) of a random variable having the distribution EC*(0i,O.5), as a function 

of 01. 
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5. Some other topics 
5.1. Alternative derivations 

As discussed in Section 2.1, the proposed model can be derived by considering a problem 
in Brownian motion. Apart from the derivation given there, here is another method to 
generate our family via this stochastic process, kindly suggested by the referee. 

Remember that, given a C- valued Brownian path without drift starting from <j) (G D), 
the point of first exit from the unit disc D has a wrapped Cauchy distribution with density 
proportional to l/\z — 0| 2 . (Although we assume <j> € D by convention, the Brownian 
motion starting at 1 / (j> has the same hitting distribution from the exterior because of the 
conformal invariance of Brownian motion.) From this it follows that another derivation 
of the proposed model (3) is given in the following theorem. 

Theorem 6. Let Z\ and Z2 be the points of first exit from the unit disc for two inde- 
pendent C-valued Brownian motions without drift starting from <p\ and 02, respectively, 
where cf>i,<f>2 £ D. Then the conditional distribution of Z\ (= Z) given Z\ = Z2 is given 
by density (3). 

Remark 1. The referee embedded the above in a generalized formulation leading to 
densities of the form 

/(z^n f-rxS ' zedD;<f> 1 ,...,<l> d eD, (10) 
j=i \ z Qj] 

many of whose properties are also readily obtainable using the residue theorem. 

In addition to the derivations using Brownian motion, there are other methods to 
generate our family. The following result shows that the proposed model appears in a 
Baycsian context. 

Theorem 7. Let Q\v have the wrapped Cauchy distribution WC(v, pi) with known p±. 
Assume that the prior distribution of v is distributed as the wrapped Cauchy WC(/i2,P2)- 
Then the posterior distribution of v given = /xi has density (2). 

The above derivation enables us to generate random samples from our model using 
the Markov chain Monte Carlo method; see Kato and Jones [10], Section 5.3, for details. 

Remark 2. In Theorem 7, the marginal distribution of Q is given by the wrapped 
Cauchy distribution, WC(/Xi + ^2, PiPz)- 

Before we discuss a third method to derive our model, here we briefly recall a known 
property of the wrapped Cauchy distribution or the Poisson kernel as follows. 
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Theorem 8 (Rudin [18], Theorem 11.9). Assume D is the closed unit disc in the 
complex plane. Suppose a function u is continuous on D and harmonic in D, and suppose 
that \u(z) \ < oo for any z S D. Then 



f u{z)±- ] l<Pl l dz = u(<t> 1 ), fceD. 
JdD 2n\z-(/) 1 \ 2 



Note that the integrand in the left-hand side of the above equation is actually the 
product of u(z) and the density of the wrapped Cauchy C*((f>i). From the fact that 

\z-<h I 2 

is a continuous real function on D and harmonic in D, it follows that 

m= " ( f 6) l~% zedD, 
2ng{(/) 1 ) 

satisfies the definition of a density. Clearly, f(z) is equal to density (3). 

The following theorem implies that our model appears as a special case of a further 
general family of distributions. The proof is straightforward and omitted. 



Theorem 9. Let f and g be probability density functions on the circle. Assume that h 
is the convolution of f and g, namely, h(r) = (/ * g)(r) = f™ n f{9)g{r — 0) dO. Then a 
function k(8) = f{0)g{r — 9)/h(r) is a probability density function on the circle. 



Density (2) can be derived on setting r = fi2 and / and g as the densities of WC(/ii, p\) 
and WC(0, P2), respectively. In this case h(r) also has the form of the wrapped Cauchy 
density WC(/xi,/5ip2) since the wrapped Cauchy distribution has the additive property. 



5.2. Mobius transformation 

The Mobius transformation is well known in complex analysis as a conformal mapping 
which projects the unit disc onto itself. It is defined as 

M(w) =a - +l3 , w£D;aedD,/3£D. 
f3w + 1 

Although this projection is usually defined on the interior of the unit disc, here we 
extend the domain of the mapping so that the boundary of the unit disc, that is, the 
unit circle, is included. It is easy to see that the unit circle is mapped onto itself via 
the transformation, namely, A4(dD) = dD. In directional statistics the transformation 
appears in some papers such as McCullagh [15], Jones [6] and Kato and Jones [9]. 
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Here we study the conformal invariance properties of our family. The distribution (3) 
has a density that is a relative invariant of weight 2 in the sense that 

f(z;<j )1 ^ 2 ) = f{M(z);M(fa),M(4>2)}\M'(z)\ 2 . 

Similarly, the extended family (10) has a density that is relatively invariant with weight 
k, in the sense that f(z;<j)i,...,(f>k) = f{M(z); M((f>i), . . . ,M(4>k)}\M'{z)\ k . The relative 
invariance with weight 1 of the wrapped Cauchy distribution can be obtained by putting 
k = l. 

The proposed family (3) is not closed under the Mobius transformation except for the 
wrapped Cauchy special cases. To say the same thing in a different way, if Z follows the 
distribution (3), then the density for M(Z) is not of the form (3) except for (f>i=Q or 
4>2 = 0. This fact can be understood clearly in the following context; as seen in Theorem 6, 
our model (3) can be derived as the conditional distribution of Z\ given Z\ = Z2, where Z\ 
and Z 2 are independent wrapped Cauchy variables. Invertibility of the transformations 
M. implies that Z\ = Z 2 if and only if M.{Z\) = M.{Z-2). However, as is known as the 
Borcl paradox (e.g., Pollard [17], Section 5.5), the conditional distribution of Z\ given 
Z\ = Z 2 is generally not the same as the conditional distribution of Z\ given M.(Z\) = 
M.{Z<i)- Consequently, the Mobius transformation of the conditional distribution of Z\ 
given Z\ = Z2, which is the Mobius transformation of model (3), is not the same as 
the conditional distribution of M.{Z\) given A4(Zi) =M.{Z<x), which is model (3). In a 
similar manner one can see that the extended model (10) with k > 2 is not closed under 
the Mobius transformation. 



6. Symmetric cases 

So far, we have mostly considered the full family of distributions with densities (2) and 
(3). In this section, we focus on the symmetric special cases of the proposed model. 
Some properties, which the general family does not have, hold for the symmetric cases. 
In Section 6.1, model (2) with p\ — p 2 or p 2 + tt is discussed. In Section 6.2, we briefly 
consider another symmetric model, namely, model (2) with p\= p 2 . 



6.1. Symmetric case I: ^ti = fi2 or n\ = ^2 + 7t 

Model (2) with fix — [i 2 or \i\ = /i 2 + n is essentially the same as the distribution with 
density 

J w = KZfx rzrzr\ ts — o la — — ; , ? 



27t(l + pip 2 ) 1 + Pi- 2pi cos(6» - p) 1 + p\ - 2/9 2 cos(<9 - p) ^ 

-7t< 9 < 71, 

where — 1 < p\ < 1, < p 2 < 1, and — n < p < n. Note the extension in the range of p: 
this density corresponds to (2) with pi = p 2 when pi > in (11), while the density is 
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equivalent to (2) with p\ = \i<x +7T when p\ < in (11). Clearly, the density is symmetric 
about 9 = p and 9 = p + 7t. 

One might, however, wish to restrict interest to the case where < p\ < 1 because 
then density (11) is unimodal. To see this, note that in this case a<j = a\ = a± = in 
(5) so that stationary points of the density occur at sin 6* = and, potentially, also at 
cos0 = —02/03 = (pi + P2XI + P1P2) I '4piP2- However, the latter expression can easily be 
proved to not be less than unity, showing that (11) is unimodal. On the other hand, the 
symmetric model (11) with negative p\ has a relationship with mixtures of two wrapped 
Cauchy distributions. The density for this submodel can be expressed as 

tins 1 1 -pj v 1 I-P2 



2n 1 + p\ - 2pi cos((9 - p) ' 2n 1 + p\ - 2p 2 cos(9 - p) ' 

where p = — pi(l — /o|)/{(l + PiP2){p2 — Pi)} - Inter alia, this representation allows random 
numbers from density (11) with p\ < to be more easily generated than they were for 
the general case in Kato and Jones [10], Section 5.3. 

6.2. Symmetric case II: pi = 

In this section we discuss the submodel of density (2) with p\ = pi = P (s [0, 1)). The 
density for this submodel reduces to 

f(a s (l-p 2 ){l + p i -2p 2 cos(p 1 -p 2 )} 
f{6) 27t(l + p2) 

x [{1 + p 2 - 2pcos(9 - pi)}{l + p 2 - 2pcos(9 - p^)}]" 1 , -n<9 <n. 

The model is symmetric about 9 = pZ and ~p + n, where ~p = arg(X^=i cos Pj + 
i Ylj=i sm Pj)- ^ s displayed in Figure 2(d), this submodel can be unimodal or bimodal de- 
pending on the choice of the parameters. The condition for unimodality can be simplified 
for this submodel as follows: 

( 2 P \ n ( 2 P 

— 2arccos < p\ — pi < 2arccos ^ 

\1 + p z J \ 1 + p z 

This result was used in Section 3.2 to find the region of Figure 1(b) in which the 
discriminant takes negative values, corresponding to unimodality, when p\ = P2 and 
pi - p 2 = 2n/3. 



7. A generalization on the sphere 

As described in Section 2.1, the proposed model (2) can be derived using Brownian mo- 
tion. By adopting a multi-dimensional Brownian motion instead of the two-dimensional 
one, we can extend model (2) to a distribution on the unit sphere. The generalized model 
is defined as follows. 
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Definition 2. Let {B t ;t > 0} be WL d -valued Brownian motion starting at Bq = p\r\\, 
where d>2, < pi < 1, m e S d and S d = {x e R d ; ||z|| =1}. Assume n = inf{t; ||B t || = 
1} and T2 = inf{i; ||i?t|| = P^ 1 }; where < p2 < 1- TTien t/ie proposed model is defined by 
the conditional distribution of B Tl given B r2 — p% r]2, where 772 G 5' d - 

For simplicity, write X = B Tl . The probability density function of this extended model 
is available, and it is given in the following theorem. 

Theorem 10. The conditional distribution of X given B T2 = p^ 1 ^ *s of the form 

1-P1P2 \\x-pim\\ d \\x-P2m\\ d y ' 

where A^-i is the surface area of S d , that is, Aj-i = 2n d ^ 2 /T(d/2) . 

See Appendix C of Kato and Jones [10] for the proof. Note that density (12) re- 
duces to the circular case (2) if d = 2, x = (cos 9, sin 6)', 771 = (cos/Ji,sin/Ji)' and 
772 = (cos/j,2, sin/^2)'- It might be appealing that density (12), which is not rotationally 
symmetric in general, can be expressed in a relatively simple form. 

When p2 = 0, the distribution becomes the so-called "exit" distribution on the sphere, 
whose density is given by 

f(x)= t~ h — ^-p , zeS d . (13) 
||ar — || rf 

We write X ~ Exitd(pi?7i) if a random variable X has density (13). This model is rotation- 
ally symmetric about x = r\\. See, for example, Durrett [3], Section 1.10, for details about 
the distribution. In particular, when d = 2, model (13) becomes the wrapped Cauchy dis- 
tribution. It is noted that this distribution is a submodel of Jones and Pewsey's [7] family 
of distributions on the sphere with density 

1 I sinhMQI'V 2 - 1 {cosh(^) + smHn^x'p} 1 ^ d 
m ^ 2V^T {d /2) Plj d ll /2 _ i{cosH ^ )} ' Xeb > () 

where P^(z) is the associated Legendre function of the first kind of degree a and order 
/3 (Gradshteyn and Ryzhik [5], Sections 8.7, 8.8). It follows from equation (8.711.1) of 
Gradshteyn and Ryzhik [5] and equation (2) of McCullagh [14] that density (14) reduces 
to density (13) if k = dlog{(l — Pi)}/2, [i = 771 and ij) = —2/d. Also, the extended 

model (12) includes the model with density 



fix) 1 ^~^ d+1 1 xeS d (15) 

— Ad-i i + p\ \\x- P imP d ' e ' (15) 



when p\= P2, 771 = 772 ■ It can be seen that this model is another submodel of Jones and 
Pewsey's [7] family (14) by putting if; = —1/d and K = dlog{(l + pi)/(l — Pi)}- Also, 
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notice that, if d = 2, density (15) corresponds to Case 2, (4), of circular submodels in 
Section 2.4. In addition to these two submodels, the multivariate model (12) contains the 
uniform distribution {p\ = p% = 0), one-point distribution (pi = and p2 1, or p% = 
and pi -4 1), and two-point distribution (pi = p2, pi ^ P2 and p\ — > 1). 

In a similar manner to Theorem 7, it is easy to prove that model (12) can be derived 
from Bayesian analysis of the exit distribution. 

Theorem 11. Let X\£ be distributed as the exit distribution Exit(pi£) with known p\. 
Suppose that the prior distribution of £ is Exit(p2?72)- Then the posterior distribution of 
£ given X = r/i is given by density (12). 

From this result and the random number generator of the exit distribution given by 
Kato [8], Section 3.3, it is possible to generate random samples from (12) through the 
Markov chain Monte Carlo method in a similar manner to Kato and Jones [10], Sec- 
tion 5.3. 
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